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Abstract 

The probability that an atom of hydrogen, initially in the f)g{m = 4) state, 
survives a 90-fs pulse of 620 nm wavelength is calculated both by direct inte- 
gration of the time-dependent Schrodinger equation and by a Floquet calcu- 
lation. The two methods give virtually identical results. The survival proba- 
bility calculated for a one-electron model of neon, for the same initial state, 
pulse duration and wavelength, is in fair quantitative agreement with the 
experimental data of van Druten et al. [Phys. Rev. A 55, 622 (1997)]. 
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As is well known, Fermi's golden rule predicts that an atom irradiated by a weak laser 
field decays by photoionization at a rate proportional to the intensity of the field if the 
photons energy exceeds the binding energy. This simple perturbative result is true, how- 
ever, only at low intensity. It is now widely accepted that the rate of ionization in an 
intense high-frequency laser field increases less rapidly with intensity than a linear power 
law, and may even decrease above a certain critical intensity (before increasing again at 
ultra-high intensity when the coupling with the magnetic field of the incident wave becomes 
important). This remarkable reduction in the photoionization rate is often referred to as 
"adiabatic stabilization", even though the atom never becomes completely stable against 
photoionization, because the quasi-stationary wave packet describing the atom at high in- 
tensity can be traced adiabatically to its field-free, initial state; in particular, its formation 
does not depend critically on fast variations of the intensity 

Adiabatic stabilization was first studied in depth in the late 1980s, by M. Gavrila, M. 
Pont and co-workers 0, and subsequently by many others. Numerical calculations focused 
initially on the stabilization of the ground state or low-excited states of hydrogen or model 
systems 0. These were extended in 1991 and later to the analysis of atoms initially in 
an excited state with a high magnetic quantum number with respect to the polarization 
direction While results interpreted by several authors as casting doubts on the existence 
of the effect were also reported 0, the numerical calculations in general showed that the 
ionization rate indeed decreases at sufficiently high intensity. In particular, the calculations 
for atoms initially prepared in an excited circular state (|m|=/~?T, — 1^1) demonstrated 
stabilization at intensities, frequencies and pulse durations that could be achieved in the 
laboratory. This prompted two experiments — the only experimental test of adiabatic 



stabilization so far — in which the photoionization of the 2p^5g{m = 4) state of neon was 
studied at a wavelength of 620 nm and intensities up to 2.3 x 10^^ W/cm^ Onset of 

stabilization was observed at high intensity, and, to the extent that a comparison could be 
made, the experimental results were found to be consistent with the results of calculations 
done at the same wavelength for the 5g{m = 4) state of hydrogen. 

The present work aims at a closer comparison of theory with the experimental data of 
Ref. The quantities measured in the experiment (the yield in photoelectrons and the 
residual populations at the end of the pulse) can be calculated, for a simplified model of the 
atom, either by direct integration of the time-dependent Schrodinger equation or by treating 
the system as being quasi-stationary and using Floquet methods. We first show that, in the 
present case, these two approaches give virtually identical results. 

Both the Floquet and the time-dependent calculations were performed in the dipole ap- 
proximation, with the field described by the vector potential A(t) = y4o(t) sincut z. For 
Aq constant and with Siegert boundary-conditions imposed on the wave function, the 
quasienergy of the (Floquet) dressed state, E, is a complex eigenvalue of a system of 
time-independent differential equations. We solved this system numerically in the veloc- 
ity gauge, using a basis of spherical harmonics and complex radial Sturmian functions P]. 
The quasienergy was also obtained by Pade-summing its Rayleigh- Schrodinger perturbative 
expansion in powers of the intensity. In the high frequency regime, the latter technique gives 
the same results as a full non-perturbative calculation for a small fraction of the cost in CPU 
time. The rate of photoionization, F = — 2Im (E/h), does not depend on the normalization 
of the Floquet wave function and is a function of the vector potential amplitude, Aq, and of 
the other parameters of the incident field. The variation of the ionization rate with intensity 
is shown in Fig. 1, for the dressed 5g{m = 3) and 5g{m = 4) states of hydrogen at 620 
nm wavelength: F increases in proportion of the intensity in weak fields, in agreement with 
the golden rule, but eventually reaches a maximum and then decreases monotonically as 
stabilization sets in. The maximum occurs at a higher intensity for m = 3 than for m = 4, 
as is discussed in Ref. [4d]. Note that in Fig. 1, F remains small compared to the inverse 
of a laser period, u/27i = 0.012 a.u., and can therefore be interpreted in the usual sense, as 
a cycle-averaged ionization rate. The process being non-resonant, one may expect that the 
atom follows the dressed states adiabatically if Aq varies slowly instead of being constant. 
Our Floquet calculations thus rehed on representing the atom by a single dressed state, i.e. 
the one that develops from the unperturbed initial state as Aq increases. Accordingly, the 
probabilities that the atom survives a (not too short) laser pulse and is left in its original 
state, Psurv) or that it is photoionized, Pion, were taken to be 



and Pion = 1 — Psurv (thus neglecting excitation of other dressed states). 

The method used for numerically integrating the time-dependent Schrodinger equation, 
for atomic hydrogen, has been described in Ref. 0. Briefly, the wave function, represented on 
a basis of spherical harmonics and (real or complex) radial Sturmian functions, is propagated 
in time using a very stable predictor-corrector algorithm based on implicit Runge-Kutta 
methods. The photoionization yield is obtained as the complement to unity of the population 
remaining bound to the nucleus at the end of the pulse, which is calculated by projection 
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on all atomic bound states with n < 25 and / < 18. When using real Sturmians, the 
basis must span a region of space large enough to contain the wave packet until the end of 
the pulse. Much smaller bases are sufficient when employing complex Sturmians, as they 
smoothly absorb the outgoing flux before it leaves the region where the atom is represented 
accurately. The results quoted below have been obtained by using complex bases, of typically 
15 spherical harmonics (with i varying from 4 to 18) and, for each partial wave, 120 complex 
radial Sturmians. We have verified that at high field intensity, 2 x 10^^ W/cm^, identical 
results (to three significant figures) are obtained with a large basis of 600 real Sturmians 
per angular momentum and 50 angular momenta. 

In order to compare the Floquet and the time-dependent methods, we calculated the 
probability that an hydrogen atom, initially in the 5g{m = 4) state with respect to the 
polarization direction, is ionized by a 120-cycle pulse. The pulse duration (90 fs, FWHM in 
intensity) was the same as in the experiment. For convenience, its envelope was taken to 
be Ao{t) = Aosm^{t/T) The results are shown in Table I. Only insignificant (< 0.0001 



%) residual populations in bound states other than the 5g state were found in the time- 
dependent calculation. Adiabatic stabilization manifests by a decrease in the ionization 
probability above about 1 x 10^'* W/cm^ peak intensity. The decrease is less pronounced in 
the ionization probability than in the ionization rate, but is still remarkable. Also noteworthy 
is the 4-digit agreement between the Floquet results and the time-dependent results. The 



agreement is better than had been previously reported for much shorter pulses ||Tl|]. It 
confirms that the one-state Floquet approximation is reliable for describing ionization in 
high frequency pulses, and makes it possible to assign, unambiguously, the reduction in the 



ionization probability observed in the time- dependent results to adiabatic stabilization |]T2 

The probability that the atom remains in the original 5g{m = 4) state is also the same 
in the two methods, since there is no net population transfer to other bound states. This 
is illustrated by Fig. 2. The figure also gives the survival probability for a pulse of same 
width as above but with a sech^ intensity profile: the probability is a little smaller for 
the sech^-pulse than for the finite-duration pulse because for the same fiuence the former 
irradiates the atom at (relatively) weak intensities for a longer time compared to the latter. 
A simple calculation shows that for a fixed pulse length and increasing fiuence, the survival 
probability tends asymptotically to a constant smaller than 1 for a sech^-pulse, assuming of 
course adiabaticity and a continuing decrease of F at high intensity. However, it decreases 
to zero for pulses with, for example, a Lorentzian profile or increases to 1 for finite-duration 
pulses. Thus adiabatic stabilization (the decrease in the ionization rate) does not necessarily 
translate into a decrease in the ionization probability. 

We now compare our theoretical results with the experimental data. We model the neon 
atom by a single electron moving in a central potential, V{r), that accounts for the partial 
screening of the nuclear attraction by the inner electronic shells. In atomic units ||T3| , 

V{r) = -i[2/o(«ir) + 2/o(«2r) + 5/i(«3r) + 1] (2) 
r 

where the parameters ai = 13.06, 0^2 = 7.2 and = 3.68 were obtained by fitting the 
energy levels of V{r) to experimental values, and 

/,M^e-"|(l-^)(f^. (3) 
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Not surprisingly, the short-range part of V{r) is of very httle importance for states with 
m = 3 or 4. The rates of ionization from the dressed 5g{m = 3) and 5g{m = 4) states of 
the model atom are within 0.02 % and 0.0005 %, respectively, from those based on a pure 
Coulombic potential (shown in Fig. 1). 

The fraction of the atomic population left in the initial state at the end of a 90-fs pulse 
is shown in Fig. 3. The surviving fraction was measured directly in the experiment, and was 
also estimated indirectly as the complement to unity of the measured ionized fraction. The 
two sets of data are consistent and do not follow the rapid decrease at high fiuence that a 
naive application of the golden rule would predict (see the dotted hne in the figure, or Fig. 
8 of Ref. H^). The solid line in Fig. 3 represents the result of the Floquet calculation for the 
sech^-pulse of above. We assumed a uniform distribution of intensity since in the experiment 
the target atoms were prepared in the initial excited state only in the central part of the 620 
nm beam. We also assumed that the system was initially in a pure 5g{m = 4) state while in 
the experiment the atoms were initially in a superposition of several correlated states with a 
small (at most 13 %) admixture of states having a 5g{m = 3) character. The importance of 
this admixture was tested by repeating the calculation for an incoherent superposition of the 
5g(m = 4) and 5g{m = 3) states in a 87 % / 13 % ratio. The resulting survival probability, 
represented by the dashed line, is smaller than for the pure 5g{m = 4) state, since ionization 
proceeds faster from the m = 3 state, but the difference is minor. As seen from the figure, 
the calculated probabilities and the data are in fair quantitative agreement. This, of course, 
tails well with the conclusion of van Druten el al. 0, that the large survival probability at 
high fiuence observed in the experiment can be attributed to adiabatic stabilization. 

The authors thank J. F. McCann for his critical reading of the manuscript. B.P. is 
supported by the "Fonds National de la Recherche Scientifique" . 
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TABLES 



TABLE I. Probability that an atom of hydrogen, initially in the 5g{m = 4) state is ionized by 

the 120-cycle laser pulse described in the text, as calculated either by full numerical integration 
of the time-dependent Schrodinger equation or by Floquet methods. The wavelength is 620 nm 
and the polarization is linear. The reduced mass of the system was taken to be 1 a.u. in both 

calculations. 

Peak Intensity Time-dependent Floquet 



0.25 X 10^^ W/cm^ 



0.10103 

0.14122 
0.16035 
0.15596 
0.14705 



0.10105 
0.14126 
0.16042 
0.15599 
0.14705 



0.50 
1.00 
1.50 
2.00 
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FIG. 1. Variation with intensity of the rate of photoionization of atomic hydrogen from the 

5g{m = 3) state (dashed Hne) or the 5g{m = 4) state (soHd Une) at 620 nm wavelength. The field 
is linearly polarized along the axis of quantization of the angular momentum. The rate is expressed 
in atomic units. 
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FIG. 2. Probability that an atom of hydrogen initially in the 5g{m = 4) state of hydrogen is in 
the same state at the end of a pulse of 620 nm wavelength and 90 fs duration (FWHM in intensity), 
as a function the pulse fluence. The field is linearly polarized along the axis of quantization of 
the angular momentum. Solid line: results of the Floquet calculation for the finite-duration pulse 
defined in the text. Squares: results of the time-dependent calculation for the finite-duration pulse. 
Dashed line: Floquet results for a pulse with a sech^-profile in intensity. 
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FIG. 3. Surviving fraction of the 2p^5g population. Open circles: surviving fraction measured 
by van Druten et al. [7]. Solid circles: complement to unity of the ionized fraction [7]. Solid curve: 
probability calculated non-perturbatively for a pure bg{m = 4) initial state. Dashed curve: same 
as the solid curve, but for an incoherent superposition of ^g{m = 4) state (87 %) and ^g{m = 3) 
state (13 %). Dotted curve: prediction of first-order perturbation theory for a pure bg{m = 4) 
initial state. 
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